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CHARM DALITZ PLOT ANALYSIS FORMALISM 
AND RESULTS (EXPANDED RPP-2004 VERSION) 

Written November 2003 by D. Asner (University of Pittsburgh) 

Introduction: Charm meson decay dynamics have been stud- 
ied extensively over the last decade. Recent studies of multi- 
body decays of charm mesons probe a variety of physics includ- 
ing doubly-Cabibbo suppressed decays [1], searches for CP vi- 
^ ; olation [1,2,3], the properties of established light mesons [4,5,6] 

q \ and the properties of nn [1,6,7] and Kn [8] S-wave states. 

^ | Future studies could improve sensitivity to D°-D° mixing [9]. 

\ Weak nonleptonic decays of charm mesons are expected 

Q . to proceed dominantly through resonant two-body decays in 

0\ , several theoretical models [10]; see Ref. [11] for a review of res- 

onance phenomenology. These amplitudes are calculated with 
the Dalitz plot analysis technique [12], which uses the mininum 



> 



number of independent observable quantities. For three-body 
O | final states when the parent particle is a scalar, the decay 

2 ! rate [13] is 

O 



\M\ 2 dm 2 12 dm 2 23 , (1) 



(2yr) 3 32v / i3 

Ph. where ma is the invariant mass of i — j and the coefficient of 

^ . the amplitude includes all kinematic factors. The scatter plot 

in m| 2 versus m\^ is called a Dalitz plot. If |.M| 2 is constant 
^ ■ the allowed region of the plot will be populated uniformly with 

events. Any variation in \M.\ over the Dalitz plot is due to 
dynamical rather than kinematical effects. 

Formalism: The amplitude of the process, D — ► rc, r — ► ab, is 
given by 

M r (L, m ab , m bc ) = ^2 ( ab \ r \) T r ( m ab) (cr\\D) (2) 

A 

= Z (L,p, q) Sf (\p\) Bl (|$|) T r (m ab ) , 

where the sum is over the helicity states A of the intermediate 
resonance particle r, a and b are the daughter particles of 
the resonance r, c is the spectator particle, L is the orbital 
angular momentum between r and c, p is the momentum of 
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c in the r rest frame and q is momentum of a in the r 
rest frame, Z describes the angular distribution of final state 
particles, and B r L are the barrier factors for the production 
of r — c and a — b, respectively, with angular momentum L 
and T r is the dynamical function describing the resonance r. 
Usually the resonances are modeled with a Breit-Wigner and 
the nonresonant contribution D — > abc is parameterized as an 
S-wave with no variation in magnitude or phase across the 
Dalitz plot. Some more recent analyses have used the X-matrix 
formalism [14] with the P-vector approximation [15] to describe 
the 7T7T S-wave. 

Barrier Factor Bj_, : The maximum angular momentum L in 
a strong decay is limited by the linear momentum q. Decay 
particles moving slowly with an impact parameter (meson 
radius) r of order 1 fm have difficulty generating sufficient 
angular momentum to conserve the spin of the resonance. The 
Blatt-Weisskopf [16,17] functions B^, given in Table 1, weight 
the reaction amplitudes to account for this spin-dependent 
effect. These functions are normalized to give Bl = 1 for 
z = (\q\ r) 2 = 1. Another common formulation B' L , also in 
Table 1, is normalized to give B' L = 1 for z = zq = (|gb|^) 2 
where qo is the value of q when m a h = m r . 



Table 1: Blatt-Weisskopf barrier factors. 



L 



B' L (q, go) 







1 



1 



1 






where z = \q\ r and zq = \qo\r 
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Angular Distribution: Z depends on the spin L of resonance 
r. The Zemach formalism [18] describes the angular distribu- 
tions in terms of Legendre polynomials but is only valid for 
reactions where a, b and c are spin-0. For final state particles 
with non-zero spin, the helicity formalism is required [19]. The 
sum over helicity states yields angular distributions that are 
proportional to Legendre polynomials when transversality is 
enforced (m 2 = m 2 ab rather than m 2 = m 2 in the helicity sum 
denominator, see Eq. (4)). 

Vector Intermediate States: The angular and barrier factors 
of Eq. (2) for vector resonances are 

(cr x \D) = [f + (P D + P c )^ + /_ (P D - P c ) J g D e^ 
(ab\r x ) = g r er (P« - A), , (3) 



where Pi is the four- momentum of particle i, e is the polarization 

4) 



vector associated with each decay vertex, /+ = f + (m 2 b ) and 



/_ = f-{fn 2 ah ) are strong interaction form factors and go and 
g r are strong interaction coupling constants. 

We evaluate the sum over helicities A = ±1, as 

S>r* = -<r+^-. (4) 



Transversality is enforced if we take the denominator of the 
second term to be m? ab rather than m 2 . This enforces a vector 
current and e^q^ = by construction and so the f-(iri^ b ) term 
does not contribute. The product of the remaining three factors 
f + g r gjj, are approximated as a constant /+(m 2 & = m 2 ) and 
g r ^D are the L — 1 Blatt-Weisskopf factors given in Table 1, 
respectively. 

The angular distribution for vector intermediate states is 
given by 

Z = (Pd + Pc), {-sT + ^/m 2 ab ) (P a - P h ) v 
= (™L - m 2 c ) + (m 2 D - m 2 c ) (m 2 - mf) /m 2 ab 
= ~2p.q, (5) 
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where the Zemach form [18], following the last equality, is only 
obtained when transversality is enforced. 

Scalar Intermediate States: There are no polarization vec- 
tors associated with the decay vertices so the process is either 
an S-wave decay or a D-wave decay. The S-wave decay has a 
uniform angular distribution and the L = Blatt-Weisskopf fac- 
tor is also constant. The D-wave decay amplitude is described 
by 



f+ (m 2 ab ) (P D+Pc ) +f_ (ml) {P D -P e ) 



9D 



T r (m ab )[g + (m 2 ab ) (P a + P b f +g _ ( m 2 ab ) (P a -P b f]g r . (6) 

The D-wave contribution is suppressed, due to small momenta, 
relative to the .S-wave and has been neglected is charm Dalitz 
plot analyses. 

Tensor Intermediate States: The angular and barrier fac- 
tors of Eq. (2) for tensor resonances are 

(cr x \D) = Ti {m 2 ab ) (P D + P c )^ (P D + P c ) v g D e^ , 
(ab\r x ) = T 2 (m 2 ab ) (P a - P b ) a (P a - P b ) p g r e*f , (7) 

where e are the polarization tensors associated with each decay 
vertex, T\{m 2 b ) and T2(m 2 b ) are strong interaction form-factors 
and go and g r are strong interaction coupling constants. 
We evaluate the sum over helicities A = ±2, ±1, [20] as 

*av ad (W^P+WPtt™) K^K"/ 3 
Z^ e A e A = 2 3 (8) 

Q Q 

where — —g^ + 9 . Transversality is enforced if we take 
the denominator of the second term to be m 2 ab rather than 
m 2 . This enforces a tensor current and q^^v — q^b 6 ^ = ® 
by construction. The product of the remaining four factors 
T\T2g r gDi are approximated as a constant Ti^rn 2 , = m 2 ) and 
g r ^D are the L = 2 Blatt-Weisskopf factors given in Table 1, 
respectively. 
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The angular distribution for tensor intermediate states is 
Z = (P D + P C )^(P D +P C ) U ± >- — 

(P a ~ P b ) a (Pa ~ P h )p 

2 2 | ( m D- m2 c) { m l- m lY ' 

™bc ~ ™ ac + ~i - 

Kb 

l/ 9 9 9 (m 2 n — m 2 ) 2 
- ml b -2m 2 D -2m 2 c + y D cJ 

\ m ac 

9 9 9 (m 2 .— m?) 2 



3(p.q) 2 -(\p\m 2 (9) 



where the Zemach form [18], following the last equality, is only 
obtained when transversality is enforced. 

Dynamical Function T r : The dynamical function T r is de- 
rived from the S'-matrix formalism. In general, the amplitude 
that a final state / couples to an initial state % through the 
unitary scattering operator S, Sfi = (f\S\i). where the scat- 
tering operator S is unitary, and satisfies SS^ = S = I. The 
transition operator T is defined by separating the probability 
that f = i yielding, 

S = I + 2iT = I + 2i{p} 1/2 f{p} 1/2 , (10) 

where / is the identity operator, T is Lorentz invariant transition 
operator, p is the diagonal phase space matrix where pa = 
2qi/m and % is the breakup momentum for decay channel i. In 
the single channel S-wave scenario S = e 2lS satisfies unitarity 
and implies 

f = -e iS sin5. (11) 
P 

There are three common formulations of the dynamical 
function. The Breit-Wigner formalism is the simplest formula- 
tion. - the first term in a Taylor expansion about a T matrix 
pole, The X-matrix formalism [14] is more general, allowing 
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more than one T matrix pole and coupled channels while 
preserving unitarity. The Flatte distribution [21] is used to 
parameterize resonances near threshold and is equivalent to a 
one-pole, two-channel X-matrix. 

Breit-Wigner Formulation: The common formulation of 
the Breit-Wigner decaying to spin-0 particles a and b is 

Tr (mat) = — 2 1 ■ p ^ (12) 

m r - m ab ~ lm r?ab (?) 

where the "mass dependent" width T a i, is 

/ \ 2L+1 / x 

T ab = T r U) f^V^o) 2 . (13) 
\%J \m ah J 

A Breit-Wigner parametrization best describes isolated, non- 
overlapping resonances far from the threshold of additional 
decay channels. Unitarity can be violated when the dynamical 
function is parameterized as the sum of two or more overlapping 
Breit-Wigners. The proximity of a threshold to the resonance 
shape distorts the line shape from a simple Breit-Wigner. This 
scenario is described by the Flatte formula and is discussed 
below. 

K -matrix Formalism: The T matrix can be described as 

f = [i-ikp^j 1 k, (14) 

where K is the Lorentz invariant X-matrix describing the 
scattering process and p is the phase space factor. 

Resonances appear as a sum of poles in the X-matrix 



m a T ai (m) m a T a j (m) 



For the special case of a single channel, single pole we obtain 

moT (m) 



K 
and 



— m z 



T = K(l-iK)- 1 = morM 



itiq — m 2 — imoT (m) 
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which is the relativistic Breit-Wigner formula. For the special 
case of a single channel, two poles we have 

m Q r Q (m) (m) 

9 9 ' 9 9 

m A a — m A nip — m z 

and in the limit that m a and ma are far apart relative to the 
widths we can approximate the T matrix as the sum of two 
Breit-Wigners, T(K a + Kp) » T(K a ) + T(Kp), 



m a T a (m) (m) 



777,2 _ m 2 _ i ma Y a (m) m 2 — m 2 — impTp (m) 



(15) 



In the case of two nearby resonance Eq. (15) is not valid and 
exceeds unity (violates unitarity). 

This formulation, which applies to s-channel production in 
two-body scattering ab — > cd, can be generalized to describe 
the production of resonances in processes, such as the decay of 
charm mesons. The key assumption here is that the two-body 
system described by the X-matrix does not interact with the 
rest of the final state [15]. The quality of this assumption 
varies with the production process and is appropriate for scat- 
tering experiments like ir~p — > 7r°7r°n, radiative decays such as 
(p,J/i(j — > 77T7T and semileptonic decays such as D — > KkIv. 
This assumption may be of limited validity for production pro- 
cesses such as pp — > nun or D — > nnn. In these scenarios the 
two-body Lorentz invariant amplitude, F, is given as 

Fi=(l- iKp} 1 Pj = fk~ x P (16) 

where P is the production vector which parameterizes the 
resonance production in the open channels. 

For the tttt S-wave the five channels relevant for D decays, 
are tttt, KK, rjrj, rj'rj' and 4tt. In this scenario a common 
formulation of the K — matrix is 

K ( A _ J V" 9j a) 9j a) sc 1 - s s c \ s- s A /2ml 



The factor is the coupling constant of the K-matrix pole 
m a to meson channel %\ the parameters f£? and Sq c describe a 



February 7, 2008 07:43 



- 8- 



smooth part of the K-matrix elements; the multiplicative factor 
— ; Z! — — suppresses a false kinematical singularity near 

{S-8A0){l-8A0) 

the 7T7T threshold - the Adler zero; and the number 1 has units 
GeV 2 . 

The production vector, with i — 1 denoting tttv, is 

where the free parameters of the Dalitz plot fit are the pro- 
duction coupling p a , and the production vector background 
parameters /fj and Sq T . All other parameters are fixed by 
scattering experiments. 

Flatte Formalism: The scenario where another channel opens 
close to the resonance position is described by the Flatte 
formulation 

T(m ab ) = ^ 2 ~( 2~ 2V 9l +92 = m r^r- (19) 

m r - m ab- % [Pm + P292) 

This situation occurs in the nn S-wave where the /o(980) is 
near the KK threshold and in the irr] channel where the ao(980) 
also lies near KK threshold. For the ao(980) + resonance the 
relevant coupling constants are g\ = and (72 = gKK- For 
the /o(980) the relevant coupling constants are g\ = g n7V and 
92 = gKK where the charged and neutral K channels are usually 
assumed to have the same coupling constant but separate phase 
space factors due to m K + 7^ m K o. 

Branching Ratios from Dalitz Plot Fits: The fit to the 

Dalitz plot distribution using either the Breit-Wigner or the 
iiT-matrix formalism factorizes into a resonant contribution to 
the amplitude Aij and a complex coefficient, aje lS ^ , where aj 
and Sj are real. The definition of a rate of a single process, 
given a set of amplitudes cij and phases Sj is the square of the 
relevant matrix element. In this spirit, the fit fraction is usually 
defined as the integral over the Dalitz plot (m a b vs m^ c ) of a 
single amplitude squared divided by the integral over the Dalitz 
plot of the square of the coherent sum of all amplitudes, 



2 

J\aje i5 iMj\ dm 2 ab dml 



Fit Fraction^- = J 1 ~ J ~ - " J 1 ~_ ' " a0 ~' ' " hc (20) 



/ \ J2k a k ei5k -M k \ dm 2 ab dm 



2 

be 
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where Mj is defined by Eq. (2) and described in Ref. [22]. The 
sum of the fit fractions for all components will in general not 
be unity due to interference. 

Reconstruction Efficiency and Resolution: The efficiency 
for reconstructing an event as a function of position on the 
Dalitz plot is in general non-uniform. Typically, a signal Monte 
Carlo sample, utilizing full GEANT [23] detector simulation 
generated with a uniform distribution in phase space is used to 
determine the efficiency. The variation in efficiency across the 
Dalitz plot varies with experiment and decay mode. 

Finite detector resolution can usually be safely neglected as 
most resonances are comparatively broad. Notable exceptions 
where detector resolution effects must be modeled are (p — > 
K + K~ and uj — > n + n~. Additionally, the momenta of a, b 
and c can recalculated with a D mass constraint. This forces 
the kinematical boundaries of the Dalitz plot to be strictly 
respected. 

Background Parametrization: The contribution of back- 
ground to the charm samples varies by experiment and final 
state. The background naturally falls into four categories: (i) 
purely combinatoric background containing no resonances, (ii) 
combinatoric background containing intermediate resonances, 
such as a real K*~ or p, plus additional random particles, (iii) 
mistagged decays such as a real D incorrectly identified as D° 
and (iv) misidentified D daughters such as D + — > tv^tv + tv + or 
D+ — > K~K + n + reconstructed as D + — > K~n + n + . 

The contribution from combinatoric background resonances 
is distinct from the resonances in the signal because the former 
do not interfere with the latter since they are not from true 
D's. The usual identification tag of the initial particle as a D° 
or a D° is the charge of the distinctive slow pion in the decay 
sequence D* + — > D°7T + or D*~ — > D 7t~. Another possibility 
is the identification of one of the D's from -0(3770) — > D°D°. 
The mistagged background is subtle and may be mistakenly 
enumerated in the signal fraction determined by a D° mass 
fit. Mistagged decays contain true D 's and so the resonances 
in the mistagged sample exhibit interference on the Dalitz plot. 
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Background from mis-identified D daughters - if present - is the 
most problematic. 

Experimental Results: The Dalitz plot analysis technique 
has been applied to the decays D — > rc, r — > ab where the 
decay products a, b and c are K or tv and the intermediate 
state r is a scalar, vector or tensor meson. More generally the 
decay products could also be the pseudo-scalar 77 or 77' mesons 
or narrow vector uj or mesons. The set of charm Dalitz plot 
analyses reported by experiments are listed in Table 2. 

Table 2: Reported Dalitz Plot Analyses. 



Decay Experiment (s) 

D° -> ivsTT+TT- Mark II [24], Mark III [25] 

E691 [27], E687 [26,29] 
ARGUS [28], CLEO[l] 
D° -> K-tt+tt Mark III [25] , E687 [29] 

E691 [27], CLEO [22] 
L>° -> K°K+n- BABAR 
D° -> K°K-n+ BABAR 
D° -> tt+tt-tt CLEO [3] 

L>° -> K S K+K- BABAR [30] 

D+ -> K-7T+7T+ Mark III [25], E687 

E691 [27], E791 [8] 
£>+ -> ZV+tt Mark III [25] 

D+ -> tt+tt+tt- E687 [4], E791 [7], FOCUS [6] 
D+ -> K+K-7V+ E687 [32] , FOCUS [33] 
L>+ -> K+K-7V+ E687 [32] , FOCUS [33] 
D+ -> tt+tt+tt- E687 [4], E791 [5], FOCUS [6] 



7r 7r — Several experiments have analyzed the 
decay D° -> K^tt+tt" ". The earliest analyses, by Mark II [24], 
Mark III [25], and E687 [26], assumed only two interme- 
diate resonances, K° s p°, iT (892)- vr+, and a significant non- 
resonant component. Additional resonances were considered by 
E691 [27] but were not found to be statistically significant. AR- 
GUS [28] and E687 [29], with more events, fit the Dalitz plot 
with six intermediate resonances: if *(892)~7r+, ifQ(1430)~7r + , 
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K° s p°, K%f (975), K%f 2 (1270), and tfg/ (1400). The nonres- 
onant contribution was negligible. The early and later E687 
results [26, 29] were consistent under similar assumptions. 
The most precise results are from CLEO [1], which includes 
three additional resonances: K° s uj,K* (1680)- tt+ and the dou- 
bly Cabibbo-suppressed iiT*(892) + 7r~. They find a much smaller 
nonresonant contribution than did the earliest experiments. 

It is not straightforward to compare or combine results 
using different descriptions of the angular distributions, barrier 
factors, resonant parametrizations, and different sets of reso- 
nances. Some of the earlier results [25-27], did not include 
barrier factors [16, 17]. Most of the earlier results [25-27, 29] 
used the Zemach formalism [18] to describe the angular shape 
of the decay pattern, while the more recent results [28, 1] use 
the helicity formalism [19]. 

The significance of the nonresonant component in the 
smaller data samples has been attributed to the presence of 
the broad scalar resonances _K"o(1430)~ and /o(1370) that were 
later observed in the larger data samples. The observation of 
a small but significant nonresonant component in the largest 
data samples suggests the presence of additional broad scalar 
resonances, the k(800) and a (500). The CLEO analysis could 
accomodate the a(500) in lieu of the nonresonant component, 
but found no evidence for the «(800). 

D° -> tt+tt-tt and D° -> K°K+K — The only signifi- 
cant contribution to the resonant substructure of D° — > 7r + 7r~7r° 
is in the p7r channels. A small nonresonant component is ob- 
served but all other nn resonances, including the a(500), 
yielded fit fractions consistent with zero. The CLEO [3] results 
for D° — > 7T + 7r~7r° are given in Table 3. 



Table 3: Dalitz fit results of D° -> n+n-n°. 



Resonance 


Amplitude Phase(°) 


Fit fraction(%) 


P + 


1. (fixed) 0. (fixed) 


76.5 ± 1.8 ±4.8 


P° 


0.56 ±0.02 ±0.07 10±3±3 


23.9 ± 1.8 ±4.6 




0.65 ± 0.03 ± 0.04 -4 ± 3 ± 4 32.3 ± 2.1 ± 2.2 


nonresonant 1.03 ± 0.17 ± 0.31 77 ± 8 ± 11 


2.7 ± 0.9 ± 1.7 
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The BABAR [30] results for D° -> K°K + K~ are given 
in Table 4. The non-^> resonant substructure in K + K~ is 
significant. Resonant constributions from ao(980)°, ao(980) + , 
and /o(980) are observed. The nonresonant and the doubly 
Cabibbo-suppressed contributions are consistent with zero. 

Table 4: Dalitz fit results of D° -> K°K+K-. 
Resonance Phase(°) Fit fraction(%) 

K°(f> 0. (fixed) 45.4 ± 1.6 ± 1.0 

K°a (980) 109 ±5 60.9 ± 7.5 ± 13.3 

K°/o(980) -161 ±14 12.2 ±3.1 ±8.6 

a (980)+K- -53 ±4 34.3 ± 3.2 ± 6.8 

a (980)-K+ -13 ±15 3.2 ±1.9 ±0.5 

nonresonant 40 ± 44 0.4 ± 0.3 ± 0.8 



Charm Dalitz plot analyses might be useful for calibrating 
tools used in B decays: specifically, to extract a from B° — > 
n+n-n°, f3 from B° -> TC°K + K~, and 7 from B ± -> DK± 
followed by D -> K°K+K- or D -> ZV+Tr" [31]. 

D+ tt+tt+tt-; a cr(500) or / (600) — The decay 
D+ -> n+n+n- has been studied by the E687 [4], E791 [7] 
and FOCUS [6] experiments. The E687 experiment consid- 
ered the p(770)°7r+, / (980)7r+, / 2 (1270)7r+, and a nonreso- 
nant component. The E791 experiment included in addition 
/o(1370)7r+ and p(1450)°7r+. Both analyses found a very large 
fraction (~ 50%) for the nonresonant contribution, perhaps 
indicating a broad scalar contribution. E791 found the non- 
resonant amplitude to be consistent with zero if a broad 
scalar resonance was included in the fit. FOCUS analyzed 
its data sample using both the Breit-Wigner formalism and 
the X-matrix formalism. The Breit-Wigner analysis included 
/o(770),/o(980),/ 2 (1270),/o(1500), <r(500), and a nonresonant 
contribution. Applying the X-matrix formalism to the S wave 
and parameterizing the p(770) and /2(1270) with the Breit- 
Wigner functions also described the FOCUS data well. 
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None of these analyses has modeled the dynamics of the 
iv + 7v + interaction. Consideration of the 1 = 2 S-wave and D- 
wave phase shifts, also measured in % + p — > Tv + Tv + n [34], could 
affect the 7r + n~ S-wave result. 

E791 finds additional evidence that the low mass im feature 
is resonant by examining the phase of the nn amplitude in the 
vicinity of the reported a(500) mass. A phase variation with 
invariant tvtv mass is consistent with a resonant contribution 
[35]. 

Table 5 gives the parameters of the a(500) determined in 
charm Dalitz plot analyses. A consistent relative phase between 
the a - (500) and p(770) resonances is observed. 

Table 5: Parameters of the ex (500) resonance. 
The amplitude and phase are relative to the 
p(770). 



Experiment E791 [7] 

Decay Mode D + — > n + n + n~ 
Amplitude 1.17 ± 0.13 ± 0.06 
Phase(°) 205.7 ±8.0 ±5.2 
m(MeV/c 2 ) 478^ ± 17 
r(MeV/c 2 ) 324+^ ± 2 l 



CLEO [1] FOCUS [6] 

D° -> K° s n+n- D+ -> n+n+n- 
0.57 ±0.13 
214 ± 11 200 ±31 

513 ±32 443 ±27 

335 ± 67 443 ± 80 



D+ -> K-ir + iT + : a k(800) ? — Indication of a broad Ktv 
scalar intermediate resonance has been reported by E791 in 
the decay D + — ► K~-k + ti + [8]. Fitting the Dalitz plot with 
Z*(892)°7r+, Zj(1430)°7r+, F2(1430)°7r+, and F*(1680)°7r+, 
plus a constant nonresonant component, E791 finds results 
consistent with earlier results from E691 and E687 with a 
nonresonant fit fraction of over 90%. Having reconstructed 
more events than the other experiments, E791 was led to 
include an extra low-mass S-wave Kw resonance to account for 
the poor fit already seen by earlier experiments: A k(800) with 
m = 797±19±43 MeV/c 2 and T = 410±43±87 MeV/c 2 much 
improved the fits. The k(800) is now the dominant resonance 
and the nonresonant fit fraction is reduced from 90.9 ± 2.6% 
to 13.0 ± 5.8 ± 4.4%. As discussed with the a(500), the K-n+ 
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S-wave result could be affected by modeling the dynamics of 
the 1 = 2 n + n + interaction. 

E791 also modeled the Kix S-wave phase variation as a 
function of Ktx mass with the -KTq(1430) resonance only and 
a nonresonant component following the parameterization of 
LASS [36]. It was necessary to relax the unitarity constraint to 
describe the E791 data [37]. The Kit S-wave phase behavior 
in this model is consistent with the model that includes the k 
resonance. 

CLEO allowed scalar K-k resonances in the fit to D° — > 
K~n + n° [22] and D° — > Kgn + n~ [1] and observed a significant 
contribution for only Kq(1430) [38]. BABAR has analyzed 
the decay D° -> K°K-n+ and D° -> ~K°K + ir~ [30]. They 
fit the former Dalitz plot with both positively charged and 
neutral ¥*(892), K* (U30), 7^(1430), ¥*(1680) and a (980) _ , 
ao(1450)~, a2(1310)~ resonances, and a nonresonant compo- 
nent. The second Dalitz plot is fit with the identical resonances 
except for the a2(1310)~. A good fit is obtained in both cases 
without including the k. 

/o(980), / (1370) and / (1500) — The proximity of the 
KK threshold requires a coupled-channel or Flatte parametriza- 
tion [21] of the /o(980) in charm Dalitz plot analyses. The 
width of the /o(980) is poorly known. E791 used a coupled- 
channel Breit-Wigner function, following the parametrization of 
Ref. [39], to describe the / (980) in D+ -> n+n+n- [5], and 
measured m r = 977 ± 3 ± 2 MeV/c 2 , = 0.09 ± 0.01 ± 0.01, 
and gxK — 0.02 ± 0.04 ± 0.03. Results similar to these are 
desirable for input to the analysis of the Df — > K + K~n + [33], 
which includes the /o(980) and ao(980). 

The quark content of the /o(1370) and /q(1500) can perhaps 
be inferred from how they populate various Dalitz plots. The 
E791 analysis of D + — > 7r + 7r + 7r _ [7] finds a contribution from 
the / (1370) but not the / (1500). The FOCUS analysis [6] 
of this decay does not find a significant contribution from the 
/ (1370). For the D+ -> n+n+n-, E687 [4] and FOCUS [6] do 
not see the /o(1370) but do see a resonance with parameters 
similar to the /o(1500), while E791 [5] observes a nn resonance 
(m = 1434 ± 18 ± 9 MeV/c 2 and T = 172 ± 32 ± 6 MeV/c 2 ) 
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that is not consistent with either meson. BABAR has found 
no evidence for either the /o(1370) or the /o(1500) in D° — > 
K°K+K- [30], while CLEO has observed the /o(1370) in D° -> 
Kg7r + 7i~ [1]. Future analyses will present a clearer picture 
only if the same resonances and model of decay amplitudes are 
applied to all Dalitz plot fits. 

Doubly Cabibbo-Suppressed Decays — There are two 
classes of multibody doubly Cabibbo-suppressed (DCS) decays 
of charm mesons. The first consists of those in which the 
DCS and corresponding Cabbibo-favored (CF) decays pop- 
ulate distinct Dalitz plots: the pairs D° - K+ 7T 7T and 
D° -> K-n+n , or D+ -> K+n+n- and D+ -> K-n+n+, 

, ™ ™ r i i , ~> ^ + "°) 

are examples. CLEO [2J has reported j^^jjQ j^- + 0) = 

(0.43+°;}J ± 0.07)%. 

The second class consists of decays where the DCS and 
CF modes populate the same Dalitz plot: for example, D° — > 
K*-n+ and D° -> K*+n- both contribute to D° -> ^tt+tt - . 
In this case, the potential for interference of DCS and CF 
amplitudes increases the sensitivity to the DCS amplitude. 
CLEO [1] has reported the relative amplitudes and phases to be 
(7.1 ± l-3±g-|±g-|)% and (189 ± 10 ± 3± J 5 )°, respectively, corre- 

sponding to B J D " - tt^pq = (0.5 ± 0.2+g f +§?)%. 
P & fi(£)° -> X*(892)-tt+) v -0.1-o.V 

CP Violation — In the limit of CP conservation, charge 
conjugate decays will have the same Dalitz plot distribution. 
The D*^ 1 tag enables the discrimination between D° and D . 
The integrated CP violation across the Dalitz plot is determined 
from 



\M\ 2 - M 



2 



dm ab dm bc / / dm ab dm bc > 

|.M| + \M\ I J 

where M. and M. are the D° and D° Dalitz plot amplitudes. 
This expression is less sensitive to CP violation than the in- 
dividual resonant submodes reported in Ref. [40]. Table 6 
reports the results for CP violation. No evidence of CP viola- 
tion has been observed. 
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Table 6: Dalitz-plot-integrated CP violation. 

Experiment Decay mode Acp{%) 

CLEO [22] D° 

CLEO [2] D° 

CLEO [40] D° 

CLEO [3] D° 



The possibility of interference between CP-conserving and 
CP-violating amplitudes provides a more sensitive probe of CP 
violation. The constraints on the square of the CP-violating 
amplitude obtained in the resonant submodes of D° — > Kgn + n~ 
range form (3.5 to 28.4) x 1(T 4 at 95% confidence level [40]. 

Guidance for Future Analyses: It is essential that a con- 
sistent formalism be adopted by all experiments that report 
Charm Dalitz plot analyses. This is necessary for the Par- 
ticle Data Group to sensibly combine results from different 
experiments. Differences in the parametrizations of the angu- 
lar distribution Z, the barrier factors Bi and the dynamical 
function T r , as well as the set of resonances r, complicate the 
comparison of results from different experiments. 

As the sizes of the data samples increase, the dynami- 
cal model which describes the intermediate resonances must 
be improved. Replacing the sum of Breit-Wigner amplitudes 
with a unitary parametrization such as the i^T-matrix is im- 
portant. Of course, other dynamical models should be studied. 
In particular, the phase shift associated with n + n + in 

D (s) - 

K~n + n + , should not be neglected and momentum- 

dependent form-factors in Eq. (3) should be considered. 

Finally, for a consistent picture of the tttt and Kn scalar 
mesons to emerge from charm Dalitz plot analyses many decays 
must be evaluated simultaneously. This will allow resonances 
decaying to different final states to be analyzed together. The 
same set of intermediate states can be used for all final states 
and input (masses, widths, couplings, spectra) from scatter- 
ing experiments like iv~p — > 7r°7r°n, radiative decays such as 
0, J/ip — > 77T7T and semileptonic decays such as D — > K-ntv can 
be consistently incorporated. 



K-jv+tv -3.1 ±8.6 
K+n-n° +9tg 

K° s n+n- -0.9±2.rt}£g 
n+n-n +1+2 ± 9 
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